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The El strength function for 15 stable and unstable Sn even-even isotopes from A — 100 till 
A = 176 are calculated using the self-consistent microscopic theory which, in addition to the standard 
(Q)RPA approach, takes into account the single-particle continuum and the phonon coupling. Our 
analysis shows two distinct regions for which the integral characteristics of both the giant and 
pygmy resonances behave rather differently. For neutron-rich nuclei, starting from ^"^^Sn, we obtain 
a giant El resonance which significantly deviates from the widely-used systematics extrapolated from 
experimental data in the valley of /3-stability. We show that the inclusion of the phonon coupling 
is necessary for a proper description of the low-energy pygmy resonances and the corresponding 
transition densities for A<132 nuclei, while in the A > 132 region the influence of the phonon 
coupling is significantly smaller. The radiative neutron capture cross sections leading to the stable 
^^*Sn and unstable ^^^Sn, ^^°Sn nuclei are calculated with both the (Q)RPA and the beyond-(Q)RPA 
strength functions and shown to be sensitive to both the predicted low-lying strength and the 
phonon coupling contribution. The comparison with the widely-used phenomenological Generalized 
Lorentzian approach shows considerable differences both for the strength function and the radiative 
neutron capture cross section. In particular, for the neutron-rich ^^''Sn, the reaction cross section is 
found to be increased by a factor larger than 20. We conclude that the present approach may provide 
a complete and coherent description of the 7-ray strength function for astrophysics applications. In 
particular, such calculations are highly recommended for a reliable estimate of the electromagnetic 
properties of exotic nuclei. 

PACS numbers: 24.10.-i, 24.60.Dr, 24.30.Cz, 21.60.Jz 
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I. INTRODUCTION 



One of the paramount and challenging goals of mod- 
ern nuclear physics is to elaborate theoretical approaches 
with not only descriptive but also predictive abilities. 
This is of particular relevance for a proper description 
of nuclei far from the valley of stability since in this 
case only limited or no information is available. Self- 
consistency between the mean field and the effective in- 
teraction is also known to be of prime importance for a 
correct exploration of the excitation properties of unsta- 
ble nuclei. Another fundamental ingredient of the model, 
known to be important even for stable nuclei, concerns 
the inclusion of more complex configurations than those 
traditionally included in the Random Phase Approxima- 
tion (RPA) or Quasi-particle RPA (QRPA). Here the 
most realistic approaches include complex configurations 
with phonons, i.e. the coupling of single-particle de- 
grees of freedom with the phonon degrees (the so-called 
phonon coupling or PC). These approaches are referred 



in the litterature as the Quasi-particle-Phonon Model fl*] , 
the (Q)RPA-HPhonon Coupling Model ((Q)RPA-(-PC) 
3 and the Extended Theory of Finite Fermi Systems 
(ETFFS) [3^. The latter is based on the Green function 
method and includes the single- particle continuum which 
is necessary for nuclei with a nucleon separation energy 
close to zero. It has been recently generalized to include 
pairing using the quasi-particle time blocking approxima- 
tion (QTBA) [|. 

These approaches have been supplemented by consid- 
ering a self-consistent mean field, see for example or 
the self-consistency between the mean field and the ef- 
fective interaction [^-Q. The latter made it possible to 
perform the calculation with one unique set of interaction 
parameters, for example the Skyrme force [^Q, instead 
of two sets of parameters used in non-self-consistent ap- 
proaches (one for the effective interaction and another for 
the mean field). These improvements - taking into ac- 
count the single-particle continuum and self-consistency- 
are of great interest, first of all, for astrophysics applica- 
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tions, but also for nuclear data evaluation. The added- 
value of such approaches lies essentially in their larger 
predictive power which provides an increased confidence 
in the calculation of structure properties for exotic nu- 
clei, especially those with a large neutron excess and / or 
a small nucleon separation energy. However, one should 
note that all the approaches developed so far are in fact 
not fully self-consistent because, as it was noted earlier in 
tti^y s^lf" consistency only at the (Q)RPA 

level and do not include more complex configurations 
into the self-consistency conditions. This is one of the 
main reasons to use some additional procedures to ex- 
clude ghost states, in particular, the spurious isoscalar 
1~ state. This is achieved through a specific fit of the 
force parameters [1, @, the use of the so-called sub- 
traction procedure in the QTBA model l^or of the so- 
called "forced consistency" method [ll|, Il3 • Neverthe- 
less, accounting for the PC and self-consistency, beyond 
any doubt, increases the quality of the microscopic nu- 
clear theory and is absolutely necessary to describe si- 
multaneously the structure of ground and excited states 
for unstable nuclei. 

The role played by the PC in the description of the gi- 
ant resonances in stable nuclei is well-known. In partic- 
ular, the PC explains approximately 50% of the observed 
width, its gross structure and sometimes some fine struc- 
ture (e.g. for the E2 isoscalar resonance in ^°*Pb [1]). 
However, the direct influence of the PC on the giant 
resonance of unstable nuclei has been much less stud- 
ied though it is expected to be as important as it is for 
stable nuclei. So far, giant resonance characteristics for 
unstable nuclei have been studied systematically within 
the (Q)RPA only (for example, see [T3l - [l6| ) and, quite 
recently, within a general approach based on sum rules 
[l7| . Note that for the reasons given below, we will dis- 
cuss in the present paper only the case of electric dipole 
resonance. 

The impact of the PC on the so-called pygmy dipole 
resonance (PDR) , which lies in the low-energy tail of the 
El giant dipole resonance (GDR) and exhausts about 
1 - 2% of the energy- weighted sum rule (EWSR) [H, [H- 
[20j . is of particular interest. First, there is no consen- 
sus at present in our understanding of some important 
questions related to the PDR (see and a recent mini- 
review 
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Second, this resonance is known to have a 
significant impact on the radiative neutron capture rate 
of astrophysical interest [l^ [H, . The importance of 
the PDR is confirmed by the simple fact that it has been 
taken into account in all modern nuclear data libraries 
(though at a phenomenological level) in addition to the 
usual GDR [20, . The question arises for exotic nuclei 
where the phenomenological approach may fail to provide 
a reliable prediction because of the specific features of the 
PDR in such nuclei and the scarce experimental data on 
which systematics is based. For this reason, a similar 
approach should be followed as for the giant resonance 
problem for unstable nuclei, i.e., as discussed above, to 
use a reliable self-consistent theory which accounts for 



PC and the single-particle continuum in addition to the 
standard (Q)RPA. 

For all these reasons, the PDR problem has re- 
cently become a subject of intensive experimental (see 
Refs. [l^ ll^[25| and references therein) and theoretical 
(see Refs. [Hla [H [13, H, [13 and references therein) 
studies. Even if the total El strength of the PDR is 
small, if located well below the neutron separation en- 
ergy, it can significantly increase the radiative neutron 
capture cross section, especially, for neutron-rich nuclei 
[l5.. 16 . 22] ■ Different measurements suggest that some 
enhancement of the El-strength could be located at low 
energies even on stable nuclei, a feature that cannot be 
described within the (Q)RPA calculations. In partic- 
ular, the large-scale QRPA calculations of [l^ predict 
PDRs which are on average 1 to 2 MeV higher in energy 
than the observed values. Many recent calculations of 
the PDR [i, [3, [U, [H [H, [li] as well as the older ones 
[sol [sij performed within the non-self-consistent Quasi- 
particle-Phonon Model confirm the need to take into ac- 
count more complex configurations than those included 
in the (Q)RPA approach, most of all the Iplh ® phonon 
or 2 quasi-particles phonon configurations. However, 
large uncertainties in the description of the PDR (in par- 
ticular, its energy and strength) remain, especially for 
unstable nuclei, and only sound microscopic models can 
shed light on its existence, as well as its relative impor- 
tance and impact on neutron capture. For example, the 
self-consistent calculations with PC [1^ [s^l have shown 
that the complex configurations give a significant contri- 
bution to the radiative neutron capture cross section for 
the unstable ^■^^Sn. 

In practice, for a proper description of the PDR, at 
least two natural physical conditions need to be fulfilled: 
first, the energy of the 1~ spurious state must be equal 
to zero; second, the theory must describe correctly the 
mean energy Eq of the El giant resonance. Only in this 
case one may expect the theory to provide a reasonable 
quantitative prediction of the PDR integral features. In 
order to satisfy these two conditions different additional 
procedures have been used. The simplest way is to ad- 
just the isovector and isoscalar effective force parame- 
ters to obtain the correct values of Eq and the spurious 
1~ level energy [51]. This is suitable for stable nuclei for 
which Eq is experimentally available. For unstable nu- 
clei, if use is made of a Skyrme force and a self-consistent 
scheme without the subtraction procedure ^4^ , it is neces- 
sary to modify some of the Skyrme parameters to obtain 
an agreement with experiment [8|. The first attempt to 
take the single-particle continuum into account exactly 
at the RPA-I-PC level using the Green function method 
was made in Ref. [33| for magic nuclei. It was shown that 
a renormalization of the SLy4 force was necessary to ob- 
tain an agreement with experiment. It is worth noting 
that this conclusion is in accordance with the studies of 
[3^ IsBj who considered this idea from a different point 
of view. 

In our previous works [1, [2l[, [2^ we realized a self- 
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consistent version of the ETFFS(QTBA) using a dis- 
cretized single-particle continuum with different kinds of 
Skyrme forces including the SLy4 one, where the velocity 
force is considered in a local approximation (sometimes 
we call it DTBA). The latter had consequences for the 
renormalization of the interaction in order to locate the 
spurious state at zero energy. On the one hand, it is of 
great interest to obtain a general information about the 
El strength function and correspondingly about the ra- 
diative neutron capture cross section for many neutron- 
rich nuclei using the well-known SLy4 forces. On the 
other hand, it is clear that the inclusion of the single- 
particle continuum along with the PC effects for non- 
magic nuclei is still a difficult problem. For these rea- 
sons, we use here our DTBA approach to calculate the 
El strength function for the long Sn isotopic chain. 

The aim of the present work is twofold. First, we calcu- 
late the PDR and GDR in the long chain of the stable and 
unstable tin isotopes using the variant of the microscopic 
self-consistent version of the ETFFS(QTBA) which, in 
addition to the (Q)RPA approach, takes into account the 
single-particle continuum (by means of a discretization 
procedure) and phonon coupling in nuclei with pairing. 
For this part we concentrate on the description of the 
integral characteristics in order to gain an insight view 
into different trends for stable and unstable nuclei and to 
compare our results with the widely-used empirical for- 
mula. Second, in order to investigate the impact of the 
PC on the radiative neutron capture cross section in sta- 
ble and unstable nuclei we calculate them both with and 
without PC, within the same scheme of calculation based 
on the SLy4 Skyrme force or a slightly modified version 
of it. Here our main attention is paid to the PDR and 
its impact on the radiative neutron capture. 



II. SELF-CONSISTENT CALCULATION OF 
THE PDR AND GDR 

A. Method 

To date there are tens of different Skyrme parametriza- 
tions serving slightly different aims and fitting some bulk 
properties of the ground state. Here we use the SLy4 
parametrization of the Skyrme force (36j which proves 
to be rather successful in describing bulk properties of 
the ground state and some excited states within the 
(Q)RPA [13. 

The ground states are calculated within the HFB ap- 
proach using the spherical code HFBRAD [sl]- The 
residual interaction for the (Q)RPA and QTBA calcu- 
lations is derived as the second derivative of the Skyrme 
functional 37]. In our calculation, several simplifications 
are performed. Namely, since up to now the QTBA ap- 
proach is designed to use the BCS-based quasi-particle 
basis, we use the HFB approach to extract the quasi- 
particle characteristics and corresponding wave func- 
tions, i.e. the occupation numbers are treated as in 



the BCS approximation. The spin-orbit residual inter- 
action is dropped. The velocity-dependent terms of the 
Skyrme force are approximated by their Landau-Migdal 
limit 'so', '40"] though some physically sounder modifi- 
cations are included. There are two kinds of velocity- 
dependent terms: the first one is oc k^(5(r — r') and the 
second one to k^^(5(r — r')k (P-wave interaction in the 
momentum space). The averaged value over the density 
of the first term gives kp/26{r — r') while that of the 
second one is zero. Such an approximation violates the 
self-consistency and one has to correct the parameters 
of the residual interaction to put the spurious center-of- 
mass state to zero. We only change here the term which is 
proportional to tikp(5(r — r') by a given factor as we take 
this term approximately. This factor is usually around 
1.0-1.25 for the Sn chain. 

In general, the ETFFS(QTBA) accounts for the single- 
particle continuum completely at the RPA level for magic 
nuclei and includes the new effect of ground state correla- 
tions caused by the PC @. However, because of the tech- 
nical difficulties connected with the pairing specificity, 
these effects are not considered in the present calcula- 
tions. The quasi-particle energy cutoff of 100 MeV is 
used. We checked that within this approach the EWSR 
is fully exhausted (for the case without the velocity- 
dependent terms) and that the use of a larger basis did 
not bring any noticeable differences. The QTBA cal- 
culations are performed with the same basis. We use 
14-16 low-lying phonons of L = 2 — 6 multipolarity and 
normal parity. They are obtained within the (Q)RPA 
with the calculated effective interaction using the same 
quasi-particle-energy cutoff. Such a consistent method 
to calculate phonons is the reason for us to use a larger 
number of phonons than in the phenomenological ETFFS 
0. In Fig. [2] we test our numerical approximation of the 
single-particle continuum discretization for the ^'^■^Sn and 
^^^Sn magic nuclei, by comparing our RPA results with 
the exact account for the continuurnby the Green func- 
tion method at the RPA level [4l|, [i^, as described in 



Ref. [33[. It turns out that both calculations are almost 



identical which confirms that the discretization proce- 
dure adopted here is quite satisfactory. 

Usually the GDR strength function is obtained from 
the experimental photoabsorption cross section which is 
fitted by a simple Lorentzian functional. The El pho- 
toabsorption cross section is related to the strength func- 
tion S'(w) as follows cr_Bi(cj) = 4.022 w Sei{oj), where the 
photon energy uj is taken in MeV, S in fm^ MeV~^ and a 
in mb. The Lorentzian fit can be used to estimate the in- 
tegral characteristics of the giant resonance [4^, i44j . The 
mean energy Eq, the resonance width F and the maxi- 
mum value of the cross section (Tq £^re extracted from the 
calculated photoabsorption cross sections under the con- 
dition that the three lowest energy-weighted moments of 
the Lorentzian and of the theoretical curve should coin- 
cide in the considered energy interval. However, in or- 
der to obtain a more complete and universal information 
about the integral characteristics it is often better to use 
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FIG. 1: (Color online) 1 strength functions within 
QRPA (dashed curves) and QTBA (solid curves) for 
™Sn,"0Sn,"2Sn,"«Sn, ^^"51^,^^*811 isotopes. 

energy moments (see below). This seems to be more ap- 
propriate for very neutron-rich nuclei where one may ex- 
pect significant deviations from a Lorentzian-like shape 
for the cross section. The same (0-30) MeV summation 
interval is used for all considered nuclei. 

The other feature of our calculation scheme is the so- 
called subtraction procedure d, 0| which should avoid a 
PC double counting from the effective interaction. This 
procedure is a direct continuation of the phenomenologi- 
cal refinement philosophy used in the first formulations of 
the ETFFS d, HI] . Because the self-consistent relativis- 
tic (Q)RPA calculations are well fitted to experimental 
data and because the subtraction procedure, in princi- 
ple, provides the correct Eq value, which is equal to the 
(Q)RPA Eq value, the authors Q obtained a fulfillment 
of the above-mentioned conditions if the energy interval 
is properly chosen, for example (10-22.5) MeV for Z=50 
nuclei or (10-25) MeV for N=50 nuclei. However, alto- 
gether with a reasonable description of the GDR's width, 
the PC gives an additional strength in the low-energy re- 
gion (see Figs lll2p . So, this low-energy contribution of 
the strength is usually neglected in such an analysis of 
the integral characteristics based on the Lorentzian fit. 



B. Results 

Here we discuss our results for the long chain of tin 
isotopes from ^°°Sn up to ^^^Sn. Figs [T] and [2] show the 
strength functions for twelve Sn isotopes. In the follow- 
ing subsections, the dipole excitations are studied and 



FIG. 2: (Color online) 1 strength functions within 

QRPA (dashed curves) and QTBA (solid curves) in 
"2Sn,i^«Sn,"2Sn,^56s^^i66gj^^ 1763^^^ p^^. 132^^ 176^^ 

isotopes, the solid green curves show calculations within the 
continuum RPA. 

the GDR and PDR analyzed along the Sn isotopic chain. 
The general idea behind such an investigation is to un- 
derstand if there are some common trends for dipole ex- 
citations in stable and unstable isotopes on both sides of 
the /3-stability valley and if they can be described within 
one unique scheme and with one unique force like SLy4. 
One can argue that such a force may not be appropri- 
ate for unstable species but this force is among the most 
suited tool and can provides us with valuable theoretical 
findings, for example, as the non-Lorentzian shape of the 
GDR for very neutron-rich nuclei such as "'^^^Sn, ^^^Sn 
(see Figl2]), as discussed below. 

1. Dipole excitations and Giant Dipole Resonances 

FigEl shows the calculated integral characteristics of 
the dipole excitations for fifteen stable and unstable Sn 
isotopes. The smearing parameter is A = 200 keV for all 
calculations. The El response integral characteristics for 
the mean energies and dispersions are calculated using 
the standard definitions 

<E>^E,,^!^, D^J^-C!^Y, (1) 
mo v mo mo 

where the energy moments m^ for the energy interval 
Ai? = Ejnax — Ejnin are calculated as follows 

mk^ dE E''S{E). (2) 
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give very similar results which possibly means that fur- 
ther experimental investigations on the ^'^^Sn GDR may 
be needed. 

Our calculations show a noticeable difference both be- 
tween (Q)RPA and ETFFS(QTBA) approaches and sta- 
ble and unstable nuclei (Figs. [T][3|). The results for in- 
tegral characteristics (Fig. [S]) clearly show the necessity 
to take the PC into account for a proper determination 
of the GDR width. For the A = 100 - 132 nuclei, the 
PC gives rise to an increase of the width by as much as 
2 MeV as compared with the (Q)RPA predictions. This 
PC effect is also important in A > 132 nuclei though to 
a lesser extent. Fig. [3] shows two distinct regions for the 
integral characteristics. The first region corresponds to 
the stable isotopes, 116 < A < 124, for which the integral 
characteristics follow the well-known phenomenological 
systematics (e.g. Eq oc A~^^^ and F oc A^^/^). The sec- 
ond region includes unstable isotopes, A > 132, for which 
these systematics fail. Similar conclusions for mean ener- 
gies can be made out of the calculations obtained within 
the RQRPAZ approach 

Here we suggest another empirical systematics which 
describes quite well the mean energy in both regions de- 
fined above: 



Eo = 78A-^/^cos^a + UA'^/^ 



sin^a. 



(3) 



FIG. 3: (Color online) Integral characteristics of the 1~ state 
for some Sn isotopes: upper panel- the GDR mean energy 
versus the atomic mass number A in the (0-30) MeV in- 
terval; lower panel- the GDR dispersion versus A for the 
same interval. RQTBA and RQRPAZ results are taken from 
Refs. 0, EO], respectively. 



First we check our approach on the stable Sn isotopes 
(^^®Sn, ^^°Sn, ^^"'Sn) and obtain, as shown in FiglH a rea- 
sonable agreement with available experimental data [20] . 
The dipole excitation for these nuclei has a well visible 
Lorentzian-like form with parameters which may slightly 
vary depending on the adopted model [23|. Though our 
results are not fitted by a Lorentzian form and only mo- 
ments are compared, it can be seen that we have a reason- 
able agreement for both the mean energy and the width 
of the El resonance. The ^'^^Sn and ^^°Sn are the only un- 
stable tin isotopes which were probed to study the GDR 
and PDR 46] . Though the SLy4 forces works quite well 
for stable nuclei there is no proof that it is good for un- 
stable ones. In Ref. [S] along with SLy4 Skyrme forces 
we probed BSk5 ^ and SkM* [Hi for ^^^Sn. In ah 
the three cases we obtained very similar mean energy 
values and widths (i.e 14.3 MeV and 2.9 MeV, respec- 
tively, for the SLy4 force), while the experimental data 
are 16.1 ± 0.7 MeV and 4.7 ± 2.1 MeV It has to be 
noted that other theoretical approaches (see Fig.|31) 0,1131 



where a — {N — Z)/A is the neutron excess and the 
factor 12A~^/-^ describes empirically the pairing gap, see 
[5l| . One can see that this is in accordance with the re- 
sults in Fig. 6 (see below). The first term is responsible 
for the collective GDR while the second one reflects the 
leftovers of the GDR, namely the non-collective particle- 
hole excitations. This formula reflects the fact that the 
mean energy depends on the superposition of these two 
kinds of excitations, while the degree of mixing is de- 
fined by the neutron excess a. The second term becomes 
important for unstable nuclei and is correlated with the 
increased low-lying El strength. As mentioned above, 
the El strength can hardly be described by a Lorentzian 
function over the whole isotopic chain. One may try also 
to separate out the GDR and PDR in the same way as 
done in Ref. [s^l using ^'^^Sn as the benchmark in the def- 
inition of the border between these two resonances. This 
procedure is somewhat artificial and anyway for unstable 
nuclei the GDR centroid energy itself does not follow the 
systematics, nor does the mean energy taken in the whole 
(0 — 30) MeV interval. Another interesting feature is that 
the dispersion (FigI3J lower panel) is minimal at neutron 
magic numbers (N=50, 82, 126) and maximal for open 
shell nuclei. Likewise for the mean energy, the dispersion 
is quite close to the lorentzianian width for stable nuclei 
only, while the contribution of the low-lying tail to the 
dispersion for the second region is increasingly important 
with the increasing of the neutron excess. 
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2. Pygmy Dipole Resonances 

The division of the dipole excitation between the GDR 
and PDR regions seems rather artificial though the na- 
ture of the vibrations is rather different: the low-energy 
neutron and proton transitional densities are vibrating 
mainly in phase, while in the GDR energy region these 
are out of phase (see for example the case of ^°*Pb in 
Rcf. [49]). But there is no general rule for defining 
the interval of pygmy dipole excitations for a given nu- 
cleus. We can only visually outline the transition re- 
gion between the GDR and PDR. We find that for the 
whole Sn chain this region is rather well described by 
the 8 to 10 MeV interval, as illustrated in Fig. ID This 
is why we consider here the (0 — 10) MeV interval for 
all the nuclei to get a general understanding about the 
low-lying dipole strength: the resulting mean energies 
< E >= Eifi and YiB{El) values are shown in Fig. [Sj 
along with the RQTBA and RQRPAZ jM resuhs as 
well as the experimental data available [33,[53|- We com- 
pare our results with these two relativistic approaches as 
they are the only available calculations of integral charac- 
teristics for both the GDR and PDR in many isotopes of 
the Sn chain and they have been obtained within a self- 
consistent scheme. It has to be noted that experimental 
data is available for ^^^Sn, ^^^Sn and ^^'^Sn isotopes up to 
the neutron separation energy (the low-lying strength is 
mostly concentrated in the (4 — 8.5) MeV interval) while 
the "pygmy region" for ^'^"Sn and ^^^Sn isotopes is not 
indicated in Ref. [52j . 

We obtain a reasonable agreement with experiment for 
the < E > values summed over the (4 — 8.5) MeV interval 
(Fig. [5]) while the integral strength is a few times larger 
than the experimental value. A similar behavior was ob- 
served in other self-consistent calculations [13], see Fig. 5. 
We find some sort of agreement with the experimental 
^^^Sn data (for our (0 — 10) MeV interval) which gives 
an integrated strength of the PDR of about (4±3)% of 
the EWSR [11], while our calculation gives 4% with the 
PC included and 2% without. At the same time the cal- 
culated and experimental GDR mean energies are rather 
different and moreover the experimental mean energy is 
out of the general trend (Fig|31). Out of our self-consistent 
calculations for the Sn chain and other non self-consistent 
calculations as in Ref. @, we conclude that the PDR is 
very model-dependent and probably the force adopted 
needs some modifications for a simultaneous description 
of both the GDR and PDR. It is rather evident also that 
further experimental investigations are needed as well. 

Figs. |3][S] demonstrate that there is a distinct differ- 
ence for characteristics of a low-lying strength for stable 
and unstable nuclei. Namely, the PC contribution to the 
Y,B{El) values, i.e. the difference between QRPA and 
QTBA predictions, is small for nuclei in the A > 132 
region, while in the A < 132 region the PC has a rather 
important impact. Moreover, for nuclei such as ^^^Sn, 
"^Sn and ^^^Sn (Fig.4, upper panel) the PDR is al- 
most completely defined by complex configurations in the 
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FIG. 4: (Color online) For each of the stable ^^*Sn and un- 
stable the upper panel shows the strength func- 
tion obtained within QRPA (dashed curves) and QTBA (solid 
curves) in the (O-lO)MeV interval and the lower pannel the 
corresponding transitional densities for protons (red curves) 
and neutrons (black curves) summed over the indicated inter- 
vals. The smearing parameter is 20 keV. 



(4 — 8.5) MeV interval. For the < E > values we also 
have a "border" a.t A = 132: in the A > 132 region there 
is almost no PC contribution and a decrease of < E >, 
while in the A < 132 region this picture is more compli- 
cated: the PC contribution "corrects" an A-dependence 
of the < E > values which would be expected within 
the QRPA approach, and in this sense the PC effect is 
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FIG. 5: (Color online) Integral characteristics of 1~ states in 
the (O-lO)MeV interval: upper panel- the mean energy ver- 
sus the atomic mass number A\ lower panel- T.B{E1) values. 
RQTBA and RQRPAZ results are taken from Ref. 0, HI. re- 
spectively. The results for Goldhaber-Teller model are based 
on the work [s^]. Experimental results for < E > and 
'SB[E1) values are taken from Refs. [30l. [s^. respectively. 



important. 

We find that the structure of the PDR excitation spec- 
trum is very specific to each nucleus. It is hardly a col- 
lective mode and can not be described by some system- 
atics like the GDR. For example, the Goldhaber-Teller 
model adapted for neutron-rich nuclei in Ref. [13] gives 
a rather smooth A-dependent behaviour for the PDR 
mean energy and its strength (Fig. [5]) which is not pre- 
dicted by our calculations and which is not really con- 
firmed by available experimental data. This is easy to 
understand from the simple Brown-Bolsterly model for 
the (Q)RPA approach. Indeed, in this model the excita- 
tion properties are determined by the structure of single- 
particle or single-quasi-particle levels. Namely, the larger 
the difference between the neighbouring levels the more 
collective the appropriate (Q)RPA l~-level. From this 
point of view, the collectivity is determined by the struc- 
ture of single- (quasi)particle levels and, therefore, the 
PDR structure is rather specific to each nucleus. More- 



FIG. 6: (Color online) Comparison of the neutron separation 
energy with the minimal particle-hole energy {Ep -\- Eh) and 
the minimal energy {Ep -\- Eh -\- w, where lj is the energy of 
the lowest phonon) for the different Sn isotopes. See text for 
more details. 



over, recently a thorough theoretical analysis on the PDR 
collectivity in ^'^^Sn was performed in Ref. ^55'| which 
shows that the collectivity is rather weak and only a few 
particle-hole configurations contribute to the PDR. The 
authors demonstrated that such contributions are force- 
dependent and cooperative but not coherent. 

In order to understand better the PC role in the PDR 
region for stable and unstable nuclei it is useful to con- 
sider the correlation between the neutron separation en- 
ergy and the beginning of the low-energy excitation spec- 
trum. In Fig. [51 we compare the neutron separation en- 
ergy with the minimal particle-hole energy {Ep-\-Eh) and 
the minimal energy [Ep + Eh + w) (where uj is the energy 
of the lowest phonon for a given isotope), which approx- 
imately determine the "beginning" of the low-lying tail 
for the (Q)RPA and QTBA models, respectively. We see 
clearly that the smaller the neutron separation energy 
the lower the energy of the first 1~ level. One can see 
also a border at A = 132, i.e a direct correlation and 
similarity in the beginning of the spectrum within both 
the (Q)RPA and QTBA for the neutron-rich isotopes. 
Due to the small neutron separation energy in the very 
neutron-rich Sn isotopes a relatively strong low-lying tail 
of the strength function arises very naturally both within 
the QRPA and QTBA (see Figs. [IE]). So, there is no no- 
ticeable difference here between these two approaches. 
In contrast, for the lighter Sn isotopes (see Figs. Ill4p a 
considerable contribution of the PC can be observed. 

To conclude our analysis of the PDR, we consider tran- 
sition densities which is now a standard way in investi- 
gating the nature of the nuclear excitations. Recently 
this analysis was performed for some tin isotopes within 
the QRPA approach [^,[5^[53] and within the relativistic 
QTBA [iil. In Fig[ll we show our self-consistent QRPA 
and QTBA results for the stable ^^^Sn isotope and the 
unstable ^■^^Sn and ^^^Sn. We obtain a rather similar 
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behaviour for ^^^Sn and ^^^Sn isotopes within the QRPA 
and QTBA approaches: below 8 MeV (^^^Sn) and 9 MeV 
(-'^^^Sn) the proton and neutron transition densities are 
in phase in both approaches, i.e. both have an isoscalar 
character and are clearly dominated by the neutron con- 
tribution at the surface. At higher energies they show an 
isovector behaviour. Globally, our results are in accor- 
dance with the QRPA results of 5, ^7} for stable nuclei. 
However, we would like to emphasize some major differ- 
ences between the QRPA and QTBA approaches as far 
as the transition densities for the stable ^^''Sn isotope 
are concerned (FigU]). Physically this corresponds to the 
fact that the PC gives a considerable contribution to the 
low-lying strength in stable nuclei, which can clearly be 
seen in Fig. 0] and in Fig. [T] for other stable nuclei. 

In summary, we find that with the inclusion of the PC 
the low-lying tail is predominantly of isoscalar nature up 
to about 8 MeV for all considered Sn isotopes while the 
(w 8 — 10) MeV interval is a transition region towards 
the isovector type of excitation which distinguishes the 
GDR. We also conclude that the inclusion of the PC is 
necessary to explain the PDR integral properties (includ- 
ing the integrated strength) in stable isotopes. Moreover 
it is mostly the PC that contributes below the neutron 
separation energy. For the A > 132 nuclei, and especially 
for unstable neutron-rich nuclei, the PC leads essentially 
to a redistribution of the PDR strength. 



III. RADIATIVE NEUTRON CAPTURE CROSS 
SECTIONS 

The presence of the PDR in neutron-rich nuclei is of 
particular interest since, if located well below the neu- 
tron separation energy, it can significantly increase the 
radiative neutron capture cross section and affect the 
nucleosynthesis of neutron-rich nuclei by the so-called r- 
process [l^ [H, [13, [E^. Similarly, the presence of ex- 
tra strength at low energy in neutron-deficient nuclei can 
be at the origin of an increase of the radiative proton 
capture or photoproton emission that take place on the 
left side of the valley of /3-stability during the so-called 
rp-process or p-process, respectively [59*]. Since such nu- 
cleosynthesis processes involve exotic nuclei that cannot 
be produced in the laboratory (at least on the neutron- 
rich side), only self-consistent calculations can provide a 
reasonable prediction of their electromagnetic excitation 
properties. The impact of our newly-derived strength 
functions on the reaction cross section are discussed be- 
low. 



A. Comparison between QTBA and QRPA 

To estimate the impact our new QTBA strength can 
have on the radiative neutron capture rate of astrophys- 
ical interest, the neutron capture cross section is calcu- 
lated using the reaction code TALYS [60|]- The strength 



function with and without the PC are included in the 
calculation of the electromagnetic de-excitation transmis- 
sion coefficients. The resulting radiative neutron capture 
cross sections calculated with the strength functions of 
Fig. [7] are shown in Fig. [8] for the three Sn isotopes. 

In Ref. [13, we studied the "3Nd(n,7)i44Nd cross 
section using the non-self-consistent ETTFS(QTBA) 
strength function. Comparison with the QRPA version 
showed that the PC inclusion increases the cross section 
by a factor 2 and improves the agreement with experi- 
ment [6l| . Very recently, the results of the radiative neu- 
tron capture cross sections calculations within the self- 
consistent relativistic QTBA were performed for the four 
tin isotopes (62] . 

Since the electromagnetic transmission coefficient cor- 
responds to an integral overlap of the de-excitation 
strength function with the nuclear level density, only 
the strength function in a restricted energy range be- 
low the neutron separation energy play an important role 
for the estimate of the radiative neutron capture rate 
[ssl . [63| . This range corresponds to 7 energies of typi- 
cally 2 MeV< <4 MeV though it may be higher in 
neutron-deficient nuclei or just before crossing a neutron 
closed shell. Therefore if located in this energy range, 
the PDR might provide quite a large contribution to the 
radiative cross section. For neutron-rich nuclei, the mean 
PDR energy is relatively low and the integrated strength 
high (Fig. [5]), s o that the PDR contribution may become 
significant [la. [T6| . 

Here we consider three different compound Sn nuclei, 
namely, the stable ^^*Sn one and unstable ^^^^Sn and^^^Sn 
isotopes. In our earlier calculations [1^ [gJ we used the 
microscopically calculated (Q)RPA and QTBA strength 
functions which have been folded with a Lorentzian in or- 
der to reproduce the expected width of the strength func- 
tion. However, such a procedure tends to smear out the 
detailed structure of the strength function that may be of 
interest in the specific energy of relevance (as discussed 
above). For this reason, we consider here the realistic 
strength functions without any Lorentzian folding. 

We consider two variants of the QTBA calculations in 
order to compare the ETFFS(QTBA) calculations with 
and without (QTBA-ws) the subtraction procedure (see 
the end of Sect. HA). In the variant QTBA-ws, the isovec- 
tor part of the calculated effective interaction strength 
fex is renormalized in order to bring (in this approach) 
the mean energy Eq (Eq. 2) to the QTBA predicted value, 
i.e to fuimi the condition Eq{QTBA) = Eq{QTBA-ws). 
It turns out that the change of the fex value is no more 
than 10% (more precisely 8, 5 and 10% for ^^^Sn, ^-^^Sn 
and ^^"Sn, respectively). We find that the difference be- 
tween the QTBA and QTBA-ws strength functions is not 
large, both of them differing significantly from the QRPA 
predictions (see Fig. [7]). The differences between QTBA 
and QTBA-ws are essentially found in the redistribution 
of the strength. This effect is however large enough to im- 
pact the radiative neutron capture cross section as shown 
in Fig. [SI In particular, FigjS] shows that the QTBA-ws 
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gives a larger cross section compared with the QTBA for 
all the three isotopes although the GDR mean energies 
are the same for both variants. The comparison vi^ith the 
corresponding RQTBA calculations for ^^^Sn(n, 7)"^'^^Sn 
cross section [62] shows, in general, similar behaviour, ex- 
cept for specific energies like those around En — 100 keV 
for which our cross section is noticeably smaller than the 



RQTBA one. These deviations may stem from the spe- 
cific structure of the strength function in the energy range 
of relevance, as well as from the use of different nuclear 
ingredients in the cross section calculation, such as nu- 
clear level densities. 

In Refs. [H, |63|, where the microscopic strength func- 
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tion is folded with a Lorentzian function, we found that 
for the stable ^^''Sn the cross sections obtained with 
QRPA and QTBA were almost identical. Comparing 
the QTBA and QRPA strength functions (FigE]) and 
cross sections (FigUl) allows us to deduce directly the 
role of the PC (without any interference of an additional 
Lorentzian smoothing). In particular, the inclusion of 
the PC increases the cross sections by a factor 2-3 in 
the QTBA-ws case (Fig. [8]). In the latter case, the cross 
section obviously follows the strength function shown in 
Fig. El i.e. the extra low-lying strength is responsible 
for an increase of the reaction cross section (and conse- 
quently of the Maxwellian-averaged reaction rate of as- 
trophysical interest) by about a factor of 3 with respect 
to the predictions based on the HFB-I-QRPA calculation 
with BSk7 Skyrme force [1^1 • More specifically, the low- 
energy El strength originating from the PC contribu- 
tion increases the cross section at 0.1 MeV (an energy 
of relevance for the r-process nucleosynthsesis) by about 
30% for i23Sn(n,7)i24Sii for the QTBA and by a factor 
of about 2 for QTBA-ws. For the "iSn(n,7)"2Sn these 
figures are even larger and for ^*^Sn(n,7)^^'^Sn we ob- 
tained a similar effect for the QTBA-ws only. Thus, in 
this section we demonstrated the noticeable sensitivity of 
the radiative neutron capture cross section with respect 
to the model (QRPA, QTBA or QTBA-ws) as well as the 
force used (SLy4, BSk7). 



B. Comparison with phenomenological models 

The Lorentzian approach has been widely used for 
practical applications, though it suffers from shortcom- 
ings of various sorts. On the one hand, the location of 
the GDR maximum energy and width remains to be pre- 
dicted from some underlying model for each nucleus. For 
many applications, these properties have often been ob- 
tained from a droplet-type model or from experimen- 
tal systematics [201 ■ shown in Fig. [31 these esti- 
mate may differ significantly from the predictions ob- 
tained from sounder microscopic models. In addition, the 
Lorentzian model tends to overestimate the El strength 
at energies below the neutron separation energy. Dif- 
ferent parametrizations or functional forms (including in 
particular an energy- and temperature-dependent width) 
have been proposed (see e.g. [13, H^) to reconcile ex- 
perimental data in the photon or radiative neutron cap- 
ture channels, but none of the proposed closed forms can 
nowadays explain the various trends observed at low en- 
ergies. Besides the Lorentzian approach cannot provide 
any predictions on the low-energy PDR, neither on its 
presence, nor on its characteristics. For this reason, it 
is of particular interest to analyze to what extent our 
predictions based on self-consistent microscopic models 
differ from those used in practical applications. 

In Fig.[Slour results are compared with those obtained 
with the phenomenological Generalized Lorentzian 
(GLO) strength function [65|]. For the stable ^^''Sn, the 



GLO cross section is rather similar to those obtained 
within the ETFFS approach, to be exact, this cross sec- 
tion curve is just between the QTBA and QTBA-ws ones, 
although the strength functions can differ at low energies 
below the neutron separation energy (Figl?]). However, 
for neutron-rich nuclei, such as ^'^^Sn and ^^"Sn, the cross 
section obtained with the GLO strength on the one hand, 
and both the QRPA and QTBA, on the other hand, dif- 
fer, especially for ^^°Sn. As shown in Sect. II, the main 
reason lies in the ^-dependence of the integral charac- 
teristics, but also in the existence of a low-lying strength 
predicted by the microscopic models. Note that the 
GLO parameters used here for ^^°Sn correspond to the 
RIPL2 recommended systematics, i.e. Eq = 14.81 MeV, 
r = 4.47 MeV and (Tq = 341.5 mb which strongly differ 
from our microscopic predictions (see Fig. [3]). Fig. [3 also 
shows the spreading of the strength function down to the 
lowest energies, i.e in the vicinity of the neutron separa- 
tion energy, while the GLO model would only provide the 
tail of the GDR strength at these energies. These com- 
parisons demonstrate the non-applicability of the empir- 
ical systematics and the necessity to make use of self- 
consistent approaches for neutron-rich nuclei. (On the 
comparisons of Ml resonances with used systematics see 
Ref.jei].) 



IV. CONCLUSION 

The electric dipole strength function has been esti- 
mated on the basis of the ETFFS(QTBA) model which 
simultaneously takes into account the (Q)RPA configura- 
tions, the more complex Iplh (g) phonon or 2 quasiparticle 
<S) phonon configurations and the single-particle contin- 
uum. For the long chain of tin isotopes, the strength 
functions have been determined within our DTBA ap- 
proach which is a discretized self-consistent version of the 
ETFFS (QTBA). The QTBA strengths have been com- 
pared with the (Q)RPA ones which allowed us to study 
the contribution of the phonon coupling along the whole 
isotopic chain. 

Our conclusions concerning the GDR and PDR prop- 
erties clearly differ depending on the nuclear region con- 
sidered, namely the A < 132 and A > 132 regions. More 
precisely, 

1. for neutron-rich A > 132 Sn isotopes, we find, both 
within QRPA and QTBA, a significant difference 
in the A-dependence of the GDR mean energy with 
respect to the standard phenomenological system- 
atics. Our Eq. [Ogives a new phenomenological sys- 
tematics; 

2. although for all considered isotopes, the PC contri- 
bution to the GDR width is very important quan- 
titatively , its contribution to the nuclei of A > 132 
region is relatively smaller than it is for the A < 132 
nuclei; 
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3. the PC contribution to the PDR integral character- 
istics < E > and Y^B{E1) summed over the (0-10) 
MeV interval is small for the neutron-rich isotopes; 

4. the transition densities in most of the low-energy 
region are mainly of isoscalar nature both within 
the QRPA and the QTBA approaches. The PC 
contribution to the transition densities also affects 
the transition densities, especially in stable nuclei. 
Globally, the isoscalar behavior is revealed on the 
energy interval considered here (1 MeV). It is not 
found that the transition density of all individual 
peaks in the low-energy region is of isoscalar nature. 
We note also that for these reasons the QRPA can- 
not explain quantitatively the isoscalar-. -isovector 
splitting of the PDR in the stable ^'^^Ce observed 
in the {a,a'^) reaction [s^], see [2l| as well. 

Such a different manifestation of the PC for nuclei with 
A < 132 and A > 132 correlates very well with the neu- 
tron separation energy. Namely, the differences are much 
smaller for neutron-rich nuclei than they are for A < 132 
nuclei. Just due to this fact the low-energy parts of the 
strength functions in neutron-rich nuclei are rather sim- 
ilar within the QRPA and QTBA. 

The radiative neutron capture cross sections for 
i24Sn,i32Sn and ^^^Sn were calculated with the QTBA 
and QRPA strength functions and shown to be sensi- 
tive to the predicted low-lying strength. Significant de- 
viations from the phenomenological GLO approach [gl] 
are also obtained for the strength functions, and con- 
sequently for the neutron capture and photoabsorption 



cross sections, for the very neutron-rich isotope ^^"Sn. 
A direct comparison between the QTBA (including the 
PC) and GLO cross sections shows that the neutron cap- 
ture cross section on very neutron-rich nuclei may be in- 
creased by 2 order of magnitude with respect to the tradi- 
tional use of phenomenological models. Our results con- 
firm the necessity to use self-consistent microscopic mod- 
els when dealing with exotic nuclei. Therefore, nuclear 
data libraries should not recommend such phenomenolog- 
ical models in this case, but rather point towards newly- 
developed microscopic large-scale calculations. 

The PC and single-particle continuum, which have 
been included, in addition to the QRPA effects, in our 
calculations, are necessary ingredients to describe the 
electric GDR and PDR properly. Nevertheless, it is nec- 
essary to use some renormalization procedures (either by 
adjusting some interaction parameters or applying the 
subtraction method) to obtain the correct value of the 
spurious state energy. In addition, the approach still 
needs to include the self-consistency at the level of com- 
plex configurations. Therefore, further developments are 
needed (see also [2l|, where some unsolved issues in the 
PDR physics are discussed). 
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